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ABSTRACT 

The statistics of peak counts in reconstructed shear maps contain information beyond the power 
spectrum, and can improve cosmological constraints from measurements of the power spectrum alone 
if systematic errors can be controlled. We study the effect of galaxy shape measurement errors on 
predicted cosmological constraints from the statistics of shear peak counts with the Large Synoptic 
Survey Telescope (LSST). We use the LSST image simulator in combination with cosmological N- 
body simulations to model realistic shear maps for different cosmological models. We include both 
galaxy shape noise and, for the first time, measurement errors on galaxy shapes. We find that the 
measurement errors considered have relatively little impact on the constraining power of shear peak 
counts for LSST. 

Keywords: Gravitational lensing: weak 



1. INTRODUCTION 

Weak gravitational lensing (WL) by large-scale cos- 
mic structures has emerged as one of the most promis- 
ing methods to constrain th e parameters of both dark 
ene rgy and dark mat t er (e. g. lAlbrecht et all (120061) : see 
also lHoekstra fc Jainl (|2008D : lMunshi et al.1 (|2008| ) for re- 
views). The COSMOS survey has provided independent 
evidence of the accelerated expansio n of the Universe 
from cosmological WL measu rements ((Schrabbac k et aLl 
120101: iSemboloni et all 1201 Over the next decade, 
the Large Synoptic Survey Telescope (LSST) and other 
large surveys covering several thousand square degrees 
will produce galaxy catalogues of unprecedented qual- 
ity0 These surveys will provide WL datasets with an 
enormous wealth of information about structure forma- 
tion, enabling not just the study of traditional two-point 
statistics like the power spectrum, but also the extraction 
of information on non-Gaussianity. The combination of 
these measurements will help substantially tighten the 
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constraints of cosmological parameters. 

Galaxy clusters are collapsed objects that provide a 
complementary probe of cosmology to the power spec- 
trum. Indeed, the cluster mass function has long been 
considered a useful probe of cosmological models. It can 
be approximated analytically, and has a strong depen- 
dence on the cosmological parameters f2 m , as and w, 
where fl m is the fractional matter density of the uni- 
verse, as is the normalization of the matter power spec- 
trum at the length scale 8/i _1 Mpc and w is the evolution 
of the equation of state of dark energy. Clusters can 
be identified and their masses measured through several 
different techniques, including x-ray observations, the SZ 
effect, galaxy cou nts and weak gravitational lensing (see 
lAllen et al.ll201ll for an overview). The measurement 
of cluster mass using WL has the advantage that it is 
independent of the luminous and dynamic properties of 
the galaxies, and is sensitive to both the baryonic and 
dark matter components. However using WL to detect 
clusters, by searching for peaks in reconstructed lens- 
ing maps, is more problematic. Shear peaks detected in 
WL surveys are often not due to single galaxy clusters, 
but to chance alignments o f structure along the line-of- 
sight (Haman a et al.l [20041) . In addition, genuine clus- 
ters that are aligned with matter underdensities along 
the line-of-sight can be missed. Cluster counts using WL 
alone t herefore tend to hav e low completeness and low 
purity (Hamana et al. 2003). 

Peaks in WL maps are a direct observable in WL sur- 
veys and can constrain cosmology, regardless of whether 
they originate from a single galaxy cluster or a random 
superposition of matter over-densities. While harder to 
predict theoretically than the cluster mass function, they 
are observationally cleaner with fewer opportunities for 
systematic errors to complicate the interpretation of the 
measurement. Contributions from filaments and other 
chance alignments encode additional information about 
the structure of matter beyond the cluster mass function, 
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making peak counts a probe of cosmological parameters 
complementary to measurements from cluster counts. 

In the past several years, there has been a signif- 
icant increase in interest in lensing peaks and other 
closely related statistics^ Most work has concentrated 
on peak counts in maps of convergence, which are easy 
to simulate but observationally harder to reconstruct 
than maps of reduced shear (see Section [2] for defini- 
tions of convergence and reduced shear). Uain et all 
(2000) studied the probability d i stribu tion function of 
the convergence and iWang et all (|2009T ) investigated its 
cumulative version, the fractional area of "hot spots" 
on convergence maps. Both statistics are similar to 
peak counts in the high- convergence limit and have been 
shown to have useful cosmology sensitivity. The frac- 
tional area statistic is also known as Vq, one of the 
three Minkowski functionals for two-dimensional thresh- 
olded fields. Minkowski functionals are related to peaks 
and had been p ropos ed as a weak-len sing statistic by 
Sato et al.l (120011) an d iGuimaraesI ([20021 ) . More recently, 



Maturi et al.l (|2010D constructed an analytical approx- 



imation to the V2 Minkowski functional, which is the 
genus statistic and also corresponds to peak counts in 
the high-threshold limit. The full set of Minkowski func- 
tionals in the context o f WL has been stu died exten- 
sively both theoretic ally dMunshi et ah 1201 If) and in ray- 
tracing simulations (|Kratochvil et al.l 12012ft . In a dif- 
ferent approac h, peak counts hay e also been studied in 
wavelet space (|Pires et al.1 120091 ) , and found to break 
the degeneracy in (crg,f2 m ) cosmological models found 
in measurements of th e power spectrum alone. 

Preliminary studies (jMarian et al.ll2009tl2010l) that de- 
fined peaks as local density maxima were based on 2D 
projections of the 3D mass distribution in low-resolution 
N-body simulations. Weak-lensing peak counts us- 
ing ra y-traced simulation s wer e su bsequently stud- 
ied by IDietrich fc Hartlapl (120 lOf) an d IKratochvil et al 



2010) and more recently inlLil f|20~TTh and Marian et al 



2012). Based on simulations with better mass resolu- 



tion, these references revealed that low-amplitude peaks 
(which typically do not correspond to single collapsed 
dark matter halos) contain more cosmological informa- 
tion than high-amplitude peaks. It should be noted that 
the range of peak heights qualifying peaks as "low" varies 
greatly between these references. 

Several other aspects of WL peak counts have al s o been 
expl ored. WL peaks we re used bv iMarian et al.l (|2011l ) 
and iMaturi et al.1 (|201lD to predict constra ints on the 
rimordial non-Gaussianity parameter /nl- lYang et al] 
20111) studied the origin of the cosmologically important 
low peaks, and found that they are ty pically caused by 
a com binat ion of 4-8 low-mass halos. IKratochvil et all 
(2012) and IMarian et all (|2012h demonstrated that cos- 
mological constraints from peaks can be ti ghtened by 
combi ning several angular smoothing scales. iPires et al.l 
(2012) compared WL peak counts directly to two other 
commonly used non-Gaussian statistics, skewness and 
kurtosis, and found the peak counts to be superior 
in information extraction from WL maps. Finally, 



9 To our knowledge, lensing peaks were first considered as 
a probe of cosmology in t he early ray-tracing simulations by 
Uain He van W acrbckc (2000), who studied the O m — dependence of 
the peak counts. 



iVanderPlas et al.l (|2012D studied the effect of masks on 
shear peak counts and showed that using Karhunen- 
Loeve analysis can mitigate biases on peak count distri- 
butions caused by masked regions, and can also reduce 
the number of noise peaks. A comprehensive study of 
the uncertainty that the presence of masks introduces 
into peak counts in real observational situations has yet 
to be performed. 

Previous work on this subject has therefore determined 
the value of peak counts in constraining cosmology, both 
alone and in combination with other lensing measure- 
ments, where peak counts can break degeneracies in cos- 
mological parameter estimation. Attention has also been 
paid to optimizing the extraction of cosm ological infor- 
matio n. Work bv lMaturi et al.l (|2010D and IMarian et al.l 
(j 2 1 If ) has concentrated on determining the optimal fil- 
ter size and shape, and how filters of different sizes can 
be combined t o increase the information extracted from 
shear maps. IDietrich fc Ha rtlap (2010) demonstrated 
the value in redshift-dependent measurement of shear 
peaks. However, to date there has been no effort to in- 
clude measurement errors in the predictions made from 
the above simulations. 

This paper introduces a framework to produce realis- 
tic galaxies that can be used to trace the shear maps 
produced in cosmological simulations, with sizes, magni- 
tudes, redshifts and signal-to-noise properties matching 
observed distributions, and measurement errors match- 
ing expected uncertainties from a ten- year LSST survey. 
We take a forward-modeling approach where we com- 
pare a dataset produced with a particular cosmology with 
datasets produced for other cosmological models. In this 
way, we can compare the expected results from the dif- 
ferent cosmologies and, eventually, determine the best-fit 
to the data. 

This paper is organized as follows. The WL formal- 
ism and aperture mass calculation are introduced in Sec- 
tion^ To calculate these shear peak statistics, we start 
with a large sui t e of N -bod y simulation s desc ribed in 
IKratochvil et al.l ()2012l ) and lYang et al.1 (j2011l ). which 
are ray-traced in order to obtain maps of the shear and 
convergence parameters for seven different cosmological 
models, covering the cosmological parameter space in 
f2 TO , erg and w. These are described in Section^ We seed 
these shear maps with source galaxies in order to obtain 
a mock dataset for each cosmological model. Realistic 
galaxies are essential to predicting realistic constraints 
on cosmology, so we must include all effects that will im- 
pact the quality of the measurements. Measurement un- 
certainties depend largely on the signal-to-noise ratio of 
the flux of the galaxy, with fainter galaxies having lower 
signal -to-noise. We therefor e use the LSST Image Simu- 
lator (Co nnolly et al.l [20101 ) input catalogues to identify 
the intrinsic properties to be used for the source galax- 
ies, such as size, magnitude and redshift. Uncertainties in 
shape measurement are determined from a large suite of 
LSST simulations, used to model the expected errors due 
to atmospheric and instrumental effects and the residual 
contributions to galaxy shape distortion after the cllip- 
ticity of the point spread function (PSF) has been in- 
terpolated to the galaxy position and deconvoluted from 
the galaxy shape. This process is described in Section |H 
and provides a mock catalogue of galaxy shapes repre- 
sentative of that which would be obtained after 10 years 
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of LSST data, for each of the seven cosmological mod- 
els. We then calculate the aperture mass over each of 
these mock catalogues, and look for peaks in the maps of 
signal-to-noise ratio (SNR) of the aperture mass statis- 
tic. The resulting peak counts are described in Section[5l 
which allow us to make the hrst realistic predictions of 
constraints on cosmological parameters from shear peak 
statistics for LSST. We describe the process by which 
we calculate the constraints on cosmology by comparing 
the different mock datasets in Section H2 and discuss our 
results in Section [7] Finally, we summarize our work in 
Section El 

2. FORMALISM 

Photons from distant galaxies are deflected by the tidal 
gravitational field of matter along the line of sight. If the 
lensed image of a galaxy is smaller than the characteristic 
scale of the lensing potential, the distortion of the galaxy 
shape can be described by a linearized lens mapping, 
given by the Jacobian 



A= (1-k) 
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where g is the reduced shear g = . The com- 
plex shear 7 = 7!+ 172 describes the distortion of 
the galaxy shape, and the convergence, k, describes the 
magnification of the galaxy image relative to its source. 
For a full derivation of these param eters, see for exam- 
ple iBartelmann &: Schneiderl (|2001l ) . 

Of course, we cannot measure the shear parameters di- 
rectly, but must estimate them from the resulting small 
distortions in observed galaxy shapes. We parametrize 
galaxy shapes by the complex ellipticity e = t\ + ie 2 , 
where the components of ellipticity are normalized mo- 
ments of the light intensity of the object hj weighted by 
a Gaussian function W(xi,x 2 ): 
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The observed galaxy ellipticity is a combination of 
the intrinsic galaxy ellipticity ei nt and reduced shear g: 

^obs — tint 

+ g. Shape noise from the intrinsic ellipticity 
of galaxies af nt = (ef nt ) is much larger than g 2 , so to 
obtain g we can average over large numbers of galaxies 
(assuming that galaxy shapes and orientations are ran- 
dom over a large enough area of the sky). In this case, 
the observed ellipticity {e bs) = (.9)- The uncertainty in 
a measurement of g, a g , is therefore a combination of the 
galaxy shape noise and measureme nt uncertainty a 2 = 

oL* + <riL„. El- Previous work (e.g . iDietrich fc Hartlad 
l201CHMaturi et al.ll20iq|Marian et al.H201lD has consid- 
ered the impact of the shape noise <Ji n t but not the mea- 
surement uncertainty a meas . 

Matter over-densities along the line of sight will cause 
the shear field, and therefore the observed shapes of 
galaxies, to be tangentially aligned around the projected 

10 We define "error" as the residual between the measured and 
true quantity, and "uncertainty" as the standard deviation of the 
differences between the measured and true quantity. 



peak of the over-density. We can use this property of 
shear fields to reconstruct the aperture mass, which is a 
weighted sum over the tangential components of galaxy 
shapes around a point . We defin e the aperture mass at 
position Qq as in [Schneider (2005|), 



M ap (6 ) = J d'OQ(8)g t {8,8 ), 

where gt is the tangential component of reduced shear 
relative to 60 defined as 



9t( 



9 ) = -(.9i cos(2<£) + g 2 sin(2<£)). 



<p is the angle with respect to the horizontal axis between 
positions 8q and 8 in the map. Note the minus sign, and 
the factor of two in the angles (necessary because shear is 
spin-2, not a vector). Q{8) is the weighting function, and 
determines the statistical properties of M ap . In practice, 
the shear field is sampled by galaxies and we measure 
the reduced shear of these galaxies. We therefore esti- 
mate the aperture mass by summing over the tangential 
components of galaxy shapes using 



1 Ns 

M ap (8 ) = ^^Q(%i,t. 



(1) 



If the weight function Q(8) follows the expected 
shear profile of a mass peak then the aperture 
mass is a matched filter for detecting mass peaks. 
We use th e spherically s ymme tric function intro- 
duced by iSchirmer et alj (|2007l ). which follows an 
NFW ()Navarro et al.|[l996D profile with exponential cut- 
offs as x —> and x 
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Here, x — 8i/8 max , where 8 max gives the radius to 
which the filter is tuned. We use a value of b.Qarcmin. 
x c is a constant, set to 0.15, which has been empir- 
ically det ermined to be a good v alue for shear peak 
counting (jHetterscheidt et al.|[2005T )f TT l The rms disper- 
sion of M ap in the case of no lensing is determined 
fro m the dispersion of the intrinsic shape noise of galax- 
ies (jBartelmann fc Schneider 200l|), 



<r{M a 



(2) 



Providing the lensing is weak within the radius of the 
aperture, <j(M ap ) will be close to the rms dispersion in 
the presence of lensing. It can therefore be used as an 
estimate of the uncertainty of the aperture mass. We can 
calculate the noise directly from the data, and look for 
peaks in the map of SNR, 



SNR((9 ) 



)9i,t 



Q 2 (0 % )gf 



(3) 



11 However, it is not yet clear whether this 
filter shape in general is the best choice for 



value or this 
the low shear 

peaks, which have been d i scovered to be cosmolo gically important 
Dietrich fc Hartlap 2010; Kratochvil et al. 2010 ) and been shown 



to be due to proje ctions of multiple clusters ( Yang et al. 201lj) since 
the publication of Hctterschcidt et al. (2005). 
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We use a pixel size of 12.2" for this map. We define 
peaks in the SNR map as all pixels in the map above a 
certain threshold having 8-connectivity (i.e. pixels which 
are connected along the sides or by the corners). While 
there are other possible definitions of peaks, this one — 
corresponding to the definition of local maxima in a pix- 
elized map — is simple and makes few assumptions about 
the underlying nature of the peaks. 

We are working with thousands of 12 deg 2 simulated 
shear and convergence maps, with each map containing 
~ 1.5 million galaxies. Calculating the aperture mass for 
all maps is a significant computational problem, which we 
address by taking advantage of the properties of graph- 
ics processing units (CPUs). The implementation of the 
aperture mass calcu l ation on the GPU is described in 
detail in [Bard et al.1 (|2012f ). By using the GPU we can 
reduce the computation time per map from several hours 
to a few minutes. 

Previous work (jYang et al.ll2011|) has determined that 
peak counts in convergence maps contain additional in- 
formation not provided by the power spectrum alone. In 
order to make a similar determination about the infor- 
mation in peak counts in reduced shear maps, we must 
also calculate the power spectrum using a simple Fourier 
transformation. We will also use this information to con- 
strain cosmological parameters, alone and in combination 
with the peak counts, as described in Section [71 

3. SIMULATIONS 

In order to predict peak counts from different cos- 
mological models we must use a large suite of N-body 
simulations representing these models, ray-traced to pro- 
duce shear maps. The large-scale structure simulations 
and shear maps we use in this analysis were created 
with the Inspector Gadg et lensing simulation pipeline 
(jKratochvil et al] lin prepI Rl on the New York Blue su- 
percomputer, which is part of the New York Center for 
Computational Sciences at Brookhaven National Labora- 
tory/Stony Brook University. In this section we describe 
the simulations and the cosmological models we chose to 
study. 

3.1. N-body Simulations 

The N-body simulatio n s are the same as those 
used in lYang et all ([20111 ). iKratochvil et all ([20121 ) and 
lYang et all (|2012l) . and consist of a series of 80 CDM TV- 
body simulations with 512 3 particles each and a box size 
of240/i~ 1 Mpc. They were run wit h a modified ver sion of 
the public N-body code Gadget-2 (jSpringerl 120051) . The 
linear matter power spectrum, which serves as input for 
the initial conditions generator N-Gen IC associated with 
Gadget-2, was created with CAMB ([Lewis et al.1 12000( 1 
for z = 0, and scaled to the starting redshift of our simu- 
lations at z = 100 according to the linear growth factor. 

The N-body simulations cover different cosmological 
models produced in multiple runs with different random 
initial conditions. A total of 50 of the runs is available 
in the fiducial cosmology, with parameters chosen to be 
{O m = 0.26,fi A = 0.74, w = -1.0, n s = 0.96, cx 8 = 
0.798, H = 0.72}. These 50 runs all used the same input 
power spectrum, but each one is a different and strictly 
independent realization. This yields a statistically robust 
set of simulations. In each of the other six cosmological 
models one parameter was varied at a time, keeping the 



Table 1 

Parameters varied in each cosmological model and weak lensing 
map set. 



WL Map Set 


C8 


w 




n A 


# of 


Identifier 










sims 


Primary 


0.798 


-1.0 


0.26 


0.74 


45 


Auxiliary 


0.798 


-1.0 


0.26 


0.74 


5 


Om23 


0.798 


-1.0 


0.23 


0.77 


5 


Om29 


0.798 


-1.0 


0.29 


0.71 


5 


wl2 


0.798 


-1.2 


0.26 


0.74 


5 


w08 


0.798 


-0.8 


0.26 


0.74 


5 


si75 


0.750 


-1.0 


0.26 


0.74 


5 


si85 


0.850 


-1.0 


0.26 


0.74 


5 



others fixed, with the following values: il m = {0.23, 0.29} 
(while Ct\ = {0.77, 0.71} such that the universe stays 
spatially flat), w = {-0.8,-1.2}, and a$ = {0.75,0.85}. 
For each of these six non-fiducial cosmological models 
5 simulations are available, where each simulation used 
a different realization of the initial conditions. Table Q] 
lists all the cosmological models with their parameters 
and number of N-body simulations. 

The shear and convergence maps, described in more 
detail in the next subsection, were generated by mixing 
simulations with different random initial conditions, and 
by randomly rotating and shifting the simulation data 
cubes. For the maps in each non-fiducial cosmology a 
mixture of all five independent N-body runs was used. In 
the fiducial cosmology, two completely independent sets 
of maps are available. One of these sets, called hereafter 
the "auxiliary" map set, was created from the five inde- 
pendent N-body runs with the same five quasi-identica£3 
initial conditions as in the non-fiducial cosmologies. The 
second map set was created by mixing lens planes from 
the remaining larger ensemble of 45 independent N-body 
runs, and will be referred to as the "primary" map set. 
This is also reflected in Table [1] 

3.2. Weak Lensing Maps 

Our pipeline uses a standard tw o-dimensional ray- 
tracin g algorithm, as described in lHamana fc Mellierl 
([2001( 1. to create the weak lensing maps from the N- 
body simul ations. Earlier w o rk usi ng similar algorithms 
includ es iSchneider et al.l (H992T) . iWam bsaans s et al.l 
(11998ft and I Jain et all (120001 ) . We refer the re ader to 
IKratochvil et al.l ([2010Ll20Tl : lYang et al.l (pOll for the 
full description of our methodology and verification of 
the accuracy of the simulations used. 

The large-scale structure from the N-body simulations 
was output as particle positions in boxes at different red- 
shifts, starting at redshift z = 2. The particles were then 
projected perpendicularly onto planes spaced 80/i _1 Mpc 
apart in a direction parallel to the central line of sight 
of the m ap. We used the triangular shaped cloud (TSC) 
scheme (|Hocknev fc Eastwood! 119981 to place the parti- 



cles on a grid on these two-dimensional density planes; 
the particle surface density was then converted into the 

12 By "quasi-identical" , we mean that the random number seeds 
to create the initial particle distributions from the power spectra 
were kept the same across all cosmological models, but the normal- 
ization of the power spectra themselves were adjusted such as to 
yield the desired erg today in every cosmology. These adjustments 
are necessary due to the difference in growth factors between the 
models. 
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gravitational potential via the Poisson equation. The al- 
gorithm then followed light rays from the observer back 
in cosmic time. The deflection angle, as well as the weak 
lensing convergence and shear were calculated at each 
plane for each light ray. These depend on the first and 
second derivatives of the potential, respectively. Between 
the planes, the light rays traveled in straight lines. 

Shear and convergence maps, spanning 12 square de- 
grees, were created for 2048 x 2048 light rays. For sim- 
plicity, we created maps assuming the source galaxies to 
be at three fixed redshifts, z s = 1,1.5,2. Each cosmolog- 
ical model is represented by 500 such 12-square-degree 
maps for convergence and shear parameters for each of 
the three source galaxy redshifts. 

4. SOURCE GALAXIES 

In this section we describe how we characterize the 
source galaxies which we use to trace the shear field. We 
wish to make our prediction for shear peak counts as 
realistic as possible, and for that it is essential that we 
make our source galaxies as realistic as possible. The 
steps we take to create the ensemble of source galaxies 
can be summarized as follows: 

• Assign a spatial position for the galaxy. 

• Assign a redshift for the galaxy. Based on redshift, 
assign the galaxy a magnitude, size and intrinsic 
shape. 

• Add reduced shear to galaxy. Re-calculate size and 
magnitude. 

• Add reduced shear error to galaxy. 

4.1. Intrinsic Properties 

As described in Section |31 we have 500 realisations of 
maps for each of 7 different cosmological models. Each 
set consists of ray-traced maps of the lensing parameters 
7i, 72 and k in three redshift bins, at z=1.0, 1.5 and 2.0. 
We consider each of the 500 realisations of one cosmo- 
logical model to be independent observations of the sky, 
and for each map we generate an independent ensemble 
of galaxies to use as tracers of the shear field in three 
dimensions. The same source galaxies are used for all 
7 cosmologies, so we are effectively observing the same 
"sky" with all cosmological models. 

We scatter the galaxies randomly across the field, en- 
suring that we have an average galaxy density of 30 
galaxies arcmin -2 , which is roughly the expected galaxy 
density usable for weak len sing analyses for an LSST 
ten- year survey in r-band t Wittman et al.l [2009). At 
this first step, we have already limited how realistic 
we can make this study: in randomly positioning the 
source galaxies, we do not take into account the shifts 
in their apparent positions due to lensing and that the 
source galaxy positions are in reality correlated with dark 
matter halos in the simulation. We decided to neglect 
these effects because it lets us shoot light-rays back- 
wards in time through the N-body simulation indiscrim- 
inately, as opposed to having to determine which light 
ray hits a fixed galaxy position. Matching galaxy den- 
sity whlijnpjjt_shearjn^s is very difficult; see for exam- 
ple iBehroozi et al.l ([20111 ) . 



One consequence of these simplifications is that 
we neglect the ma gnification bias present in lensing 
(|Turner et al.l (fl984[ )V The magnification bias arises 
from two competing effects: i) high-shear regions mag- 
nify galaxies, thus making fainter galaxies visible in a 
flux-limited survey and adding source galaxies in those 
regions of the sky, ii) the magnification also spreads apart 
the apparent positions of the source galaxies, thus dilut- 
ing the number of galaxies in these high-shear regions. 

We anticipate that this variation in density will have 
a small impact on lensing peak counts, or at least on 
the cosmological constraints coming from lensing peak 
counts. This is because the constraints have been 
shown to be dominated by the numerous low peaks 
(|Kratochvil et al.ll2010t [Dietrich fc Hartladl2010D . which 
are to be found in low shear regions, while the magnifica- 
tion bias is most noticeable in regions of high shear. For 
the high significance lensing peaks, the primary effect of 
magnification on galaxies will be a shift in the apparent 
position of galaxies which can also shift the position of a 
peak (particularly if it is not the central peak of a clus- 
ter). Since peak counts (measured using one smoothing 
scale as done in this paper) do not measure angular cor- 
relations, a shift in position will not affect the results. 
For the central peak of a cluster, however, the dilution of 
source galaxies associated with magnification will mostly 
cause an apparent broadening of the peak, which will 
make the peak appear larger than in our simplification. 
We do not expect this to be a significant effect for the 
cosmological constraints, because high central peaks are 
by far outnumbered by the others, but the importance 
of this effect should be studied in future work. 

Next we assign each galaxy a redshift, size and mag- 
nitude, taken from a distribution obtained from the in- 
put catalogues of the LSST Image Simulator. Galax- 
ies in these catalogues have properties based on those 
produced by the Millennium dark matter simulations. 
The galaxy catalogue is complete out to an r magni- 
tude of 28, which is approximately one magnitude deeper 
than the expected depth of the full LSST ten-year sur- 
vey. These quantities have been anchored to observations 
from a compilation of d eep survey datJ^I. the DEEP2 
survey ([Coil et al.l 120041) , and data from the publicly 
available Hubble Deep Field catalogue^. A redshift 
is assigned at random to the galaxies, shown in Fig- 
ure [TJ where the dashed line represents a simple model of 
the fo rm n(z) cx z 2 e~ 2z , as described in IWittman et al.l 
(2000) and previously found to be a good fit to DEEP2 
survey data ([Coil et al.ll2004[ ). A redshift-dependent size 
and magnitude is assigned for each galaxy from the simu- 
lated input catalogues, where we define size as the prod- 
uct of the RMS of the semi-major and semi-minor axes 
of the galaxy. 

An intrinsic ellipticity is then assigned to each 
galaxy. We base the assigned el lipticity on measure- 
ments made usin g COSMOS data (jJoachimi et al.ll2f)T2t 
iLeauthaud et a l. 2007), where the intrinsic galaxy shape 
noise was found empirically to be 0.23 per reduced shear 
component. A small dependence on galaxy magnitude 
was identified in the COSMOS data, but no significant 
dependence on size or redshift was found. We next as- 

13 http:/ /astro. dur.ac.uk/~nm/pubhtml/counts/counts. html 

14 http:/ /www. stsci.edu/ftp/science/hdf/archive/v2. html 
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Figure 1. Number of galaxies as a function of redshift for the 
source galaxies. The blue dashed line shows a fit to the simple 
model n(z) oc z 2 e~ 2z . 

sign a reduced shear to the galaxy, taken from the input 
shear and convergence maps described in Section [3] and 
extrapolated to the galaxy position in RA, DEC and red- 
shift (where we linearly extrapolate the shear and con- 
vergence parameters between the three redshift planes at 
z=[1.0,1.5,2.0]). The galaxy size, magnitude and elliptic- 
ity are re-calculated to account for the effect of shearing 
and magnification. 

4.2. Measurement Error 

To assign a measurement error for the galaxy shape is 
somewhat complicated. Our aim is to obtain an error 
consistent with a measurement of galaxy shapes from a 
ten-year stack of LSST observations, using only observa- 
tions made with good seeing in the r-band. If we require 
the median seeing to be 0.66", which is an acceptabl e 
quality for weak lensing analyses (jWittman et alj |2009). 
we expect a ten-year stack to consist of approx. 100 ob- 
servations. 

We use the LSST Photon Simulator (PhoSim) to sim- 
ulate a ten-year stack of LSST observations of an area on 
the sky the size of a LSST chip (i the area of an LSST 
CCD, with approx. 13.6 arcmin 2 field of view), at differ- 
ent positions on the LSST focal plane in order to sam- 
ple the PSF as it varies across the focal plane. PhoSim 
is a high-fidelity, ray-traced end-to-end simulator of the 
LSST system. A detailed description of the syste m can 
be found inlPeterson et all (120091). [ Connolly et all (12010ft 
and IChang et al.l (|2012b[ ). Recent work bv iChang et al.l 
(2012a]) has studied in depth the impact on galaxy shape 
measurement made by the distortions introduced by the 
atmosphere and the LSST telescope itself. We wish to 
isolate the impacts of measurement and algorithmic ef- 
fects, and to evaluate the impact of these errors sepa- 
rately from the error due to galaxy intrinsic shape noise. 

Since we have already accounted for shape noise in a 
previous step in the pipeline, we use in these simula- 
tions an input catalogue consisting of circular galaxies 
with a Gaussian profile to remove any effects of shape 
noise. The magnitudes, redshifts, SNRs and spectral en- 
ergy distributions of these Gaussian galaxies are the same 
as the fully realistic galaxy distribution of the ImSim in- 
put catalogues. The advantage to this approach is that 
we can easily evaluate the measurement error without 
attempting to remove galaxy shape noise. We are not 
performing a redshift-dependent measurement, and we 
ignore for this work the potentially substantial errors in 



redshift measurement. 

We take 500 values of [71 , 72, «] at random from one 
of the simulated shear maps described in Section |3l and 
produce 500 different sheared ImSim input catalogues 
by applying a single reduced shear value [gi, in, 92, in] to 
all galaxies in an existing catalogue. For each of these 
sheared input catalogues, we produce 100 simulated im- 
ages of the same area of sky, each time with different 
atmospheric condi tions specified b y the LSST Opera- 
tions Simulator (jPinto et al.l 120091 ) selected such that 
the median seeing is 0.66". We process the resulting 
ima ges using the SourceExtr actor object-detection pack- 
age jBernn[&3rnouti[T996) . For each exposure, we use 
the stars in the field to reconstruct the PSF, which is 
interpolated to the galaxy locations using a third order 
polynomial interpolation function. The measured galaxy 
shapes are correcte d for distortions du e to the PSF using 
the popular KSB (Kai ser et al.l 119951 ) algorithm imple- 
mented in the IMCA10 pipeline. We use the KSB al- 
gorithm because it is well known in the community, and 
its strengths and weaknesses are well understood. For 
example, it is known that the process that converts cllip- 
ticity to reduced shear should be calibrated using simu- 
lations. We apply a "perfect" calibration, by applying a 
calibration factor that shifts the mean measured reduced 
shear in each of our simulated exposures to the input 
reduced shear value. The measured, PSF-corrected, cal- 
ibrated shape for each galaxy is then averaged over the 
100 atmospheric realisations, giving us an estimate of the 
stacked galaxy shape measurement. More sophisticated 
algorithms are expected to give a smaller uncerta i nty on 
galaxy shape measurement (see IKitching et al.l ()2012D 
for a summary of the performance of many current shape 
measurement methods). Despite applying a perfect cal- 
ibration, for this reason we consider the uncertainty ob- 
tained from our KSB implementation to be conservative 
for LSST. 

We compare the measured galaxy shapes to the input 
reduced shear values, and the difference between input 
and output gives the uncertainty on the reduced shear 
measurement. There is a dependence of the measurement 
uncertainty with magnitude, with fainter galaxies having 
larger uncertainties, as shown in Figure [2] This is in ac- 
cordance with the dependence of reduced shear measure- 
ment uncertainty with obje ct magnitude found in data 
from the COSMOS survey (ILeauthaud et al.ll200l . We 
account for this dependence as we assign measurement 
errors drawn from this distribution, which are added to 
the galaxy shape noise. Since we assign the noise to 
galaxies randomly, we do not consider any spatial corre- 
lations the noise may have across the field. IChang et al.l 
(|2012a|) has shown that, for a ten-year stack of LSST 
images, the spatial correlation of measurement error (in- 
cluding an imperfectly modeled PSF) is at a level com- 
parable to the statistical error on the weak lensing cor- 
relation function, around 10~ 7 . We therefore consider 
that the spatial correlations of measurement error will 
be similarly small for aperture mass statistics, and ne- 
glect it in this work. However, future work is planned 
to specifically quantify the impact of correlated error on 
peak counts. 

We also wish to investigate the dependence of measure- 

15 http:/ /www. ifa.hawaii.edu/~kaiser/imcat/ 
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merit error with galaxy shape. To do this, we therefore 
made a set of simulations identical to those described 
above, but using elliptical galaxies. However, we are 
unable to separate the intrinsic shape noise from the 
measurement error in these simulations, so we are lim- 
ited to examining the dependence of the total error on 
the measured reduced shear with the measured galaxy 
shape. Figure [3] shows the uncertainty distribution for 
input reduced shear values [ffim,<72,m] for different val- 
ues of measured galaxy shape [<7i, me0 s, <?2,meas]- There is 
no dependence with galaxy shape, and the distribution 
is remarkably flat. We do see a significant dependence 
of the uncertainty on [gi } i n , 92,in] with the raw elliptic- 
ity measurement of the galaxy, but the process of PSF 
deconvolution using KSB, and the calibration procedure 
removes this dependence. We therefore do not apply a 
shape-dependent measurement error. 



but is not negligible. We shall see in the next sections 
what impact this has on the peak counts and cosmolog- 
ical constraints. 




23 24 
magnitude 

Figure 2. Standard deviation of the errors for measured cosmic 
shear parameters gn n and (?2,in, for circular galaxies of different r 
magnitudes. The standard deviation is taken from fitting a Gaus- 
sian to the distribution of galaxy shape measurement errors for 
each magnitude bin. Black squares represent gl, blue triangles g2. 




Figure 4. Distribution of galaxy shapes \g\, for intrinsic ellipticity 
only (black dotted line) and including measurement error (red solid 
line). 



5. PEAK COUNTS 

We consider these mock galaxy shape measurements to 
be a representative sample of an LSST ten-year survey. 
We use these simulated datasets to perform the aperture 
mass calculation given in Section [21 using t he GPU im- 
plementation described in iBard et all (|2012f ) , and obtain 
peak counts for the seven cosmological models described 
in Table [T] The aperture mass is calculated using source 
galaxies with shape noise alone, and using source galax- 
ies with both shape noise and measurement error. Figure 
[5] shows the distributions of peak counts for the different 
cosmological models, where each model is sampled by 
the same galaxies (including shape noise and measure- 
ment errors) scaled to the full-sky LSST survey size. 
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Figure 3. Standard deviation of the errors for input reduced shear 
parameters gi : j n and 92, in, for elliptical galaxies of different mea- 
sured reduced shear g\ t meas and g2,meas- The standard deviation 
is taken from fitting a Gaussian to the distribution of galaxy shape 
measurement errors for each measured g bin. Black squares repre- 
sent gi, blue triangles g%. 

Figure H] shows the distribution of the assigned values 
of the reduced shear |g| for all galaxies in our sample, 
comparing the intrinsic ellipticity alone to the combina- 
tion of intrinsic ellipticity and measurement error. The 
measurement error has a much smaller contribution to 
the total galaxy shape error than the intrinsic ellipticity, 




-2-101234 

SNR threshold 

Figure 5. Peak counts above SNR threshold for different cos- 
mological models, for the aperture mass calculated using realistic 
galaxies with both intrinsic shape noise and measurement error. 
See Table [l] for details of cosmological model parameters. 

To evaluate the impact of measurement error, com- 
pared to shape noise alone, we calculate the difference be- 
tween the peak counts for the two cases. This is shown in 
Figure [H where we plot the fractional difference between 
the peak counts for intrinsic shape noise alone, compared 
to shape noise and measurement error. The difference is 
largest at very low and high peak significance, where it 
reaches up to 25%. As we might expect, the difference 
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is identical for all cosmological simulations, showing that 
measurement error should not bias constraints on cosmo- 
logical parameters in favour of one model over another. 
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SNR threshold 

Figure 6. Fractional difference between peak counts with shape 
noise alone, and peak counts with shape noise and measurement 
error, given as the fractional difference from shape noise only. Solid 
curves are the average over the 500 different maps; error bars are 
the standard deviation of the 500 maps, shown for the auxiliary 
fiducial model to indicate the level of statistical error. See Table [T] 
for details of cosmological model parameters. 
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SNR threshold 

Figure 7. Difference between peak counts in fiducial cosmology, 
and peak counts in other cosmologies for aperture mass calculated 
using realistic galaxies with both shape noise and measurement 
error, given as the fractional difference from fiducial cosmology 
counts above SNR threshold. Solid curves show the mean differ- 
ence for the 500 maps used in the measurement; error bars are 
the standard deviation, shown for the auxiliary fiducial model to 
indicate the level of statistical error. See Table [l] for details of 
cosmological model parameters. 

It is hard to distinguish between the different cosmo- 
logical models by eye in Figure [SJ but if we plot the dif- 
ference of the peak counts from the fiducial cosmology, 
as in Figure [Jj we obtain a clearer view of the character- 
istics of each cosmology. The red curve in these figures 
represents the peak counts obtained from the auxiliary 
WL map set of the fiducial model, and acts as a control 
test to be compared with the primary map set of the fidu- 
cial model belonging to the same cosmology. In Figure [7J 
it is clear that the difference between the primary and 
auxiliary map sets of the fiducial model are consistent 
within the statistical error (shown by the error bars on 
the auxiliary fiducial model) . Several of the cosmological 
models have very similar peak count distributions. In ac- 
cordance with expectations of the (as, f2 m )-degeneracy, 
the models Om23 and si75, and Om29 and si85, have very 



similar peak count profiles, which will result in our pre- 
dicted cosmological constraints exhibiting a degeneracy 
in the corresponding direction in the parameter space. 

We can also compare these peak counts to the counts 
obtained from calculating the aperture mass directly 
from the maps of reduced shear, without sampling the 
maps with source galaxies. This is the "perfect" case 
where we have perfect knowledge of the shear, and no 
noise is introduced by galaxy shape noise or measurement 
errors. It is therefore an impossible ideal, but serves as 
a useful comparison to examine how real measurements 
are affected by error. Since there is no noise in this mea- 
surement, constructing a SNR map is meaningless, and 
instead we count peaks in the map of aperture mass. The 
two quantities can be related by 

M ap {6 ) = SNR(6 )a Map 

where M ap , SNR and <xm qP are defined in Equations [1] 
[3] and respectively. Figures H] and O show the peak 
counts above aperture mass thresholds for different cos- 
mologies, and the difference of the seven cosmological 
models compared to the primary map set of the fiducial 
model, respectively. 

The addition of shape noise and measurement error 
has a significant impact on the shape of the peak counts, 
visible in a comparison of Figures [Jj and [9] For the re- 
alistic case with errors included, the overall shape of the 
peak counts will be valuable in constraining cosmological 
parameters, since the deviation from the fiducial cosmol- 
ogy is visible at all SNR levels. The addition of noise 
has impacted both the significance of the peaks, as one 
would expect, and also the shape of the peak counts of 
the different cosmologies. 
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Figure 8. Peak counts above SNR threshold for different cosmo- 
logical models for the "perfect" case where the aperture mass is 
calculated directly from the reduced shear maps. See Table [T] for 
details of cosmological model parameters. 



5.1. The Significance of Low Peaks 

In this section we discuss further the signifi cance 
of low peaks. Both Dietrich fc Hartlad (|2010f ) and 
IKratochvil et al.l ()2010 ) discovered that low peaks 
contribute the most to the cosmological constraints, 
which was later c o nfirme d and elaborated up on b y 
others lYane 



et~aT1 (I20T1: IKratochvil et al.1 (I20II: 



Maria n et al.l (|2012[ ) ; iMaturi et al.l (|2011l ); lYang et al.l 
(|2012f) . What has gone mostly unnoticed is that the 
different papers have completely different definitions 
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aperture mass threshold 

Figure 9. Difference between peak counts in fiducial cosmology, 
and peak counts in other cosmologies for the "perfect" case where 
the aperture mass is calculated directly from the reduced shear 
maps, given in % difference from fiducial cosmology counts above 
SNR threshold. Solid curves show the mean difference for the 500 
maps used in the measurement; error bars are the standard devi- 
ation, shown for the auxiliary fiducial model to indicate the level 
of statistical error. See Table [T] for details of cosmological model 
parameters. 

of the word 'low' in this context, and different pa- 
pers actually refer to completely disjoint peak ranges. 
iKratochvil et all (120101) a nd the group's follow-up works 



( Yang et all 12011 
Group A, define 



2012), subsequently referred to as 
low peaks as having an SNR be- 
depending on the publication. 



tween 0-2 or 1-3.5 

IDietrich fc Hartfapl (|2010|) define low peaks as lying in 
the range 3.25-4.5<7, such that their entire range is higher 
than most of the previous gr oup's papers. The h ierarchi- 
cal peak finding algorithm of Ma rian et al.l (|2012f) is over- 
whelmed by the number of peaks below SNR^ 3(7 and 
breaks down, so these authors also restrict themselves to 
a range abov e SNR ^ 3cr in their peak detection, while 
IMaturi eT al. (2011) conclude that to constrain the pa- 
rameter f^L of primordial non-Gaussianity, only peaks 
with SNR> 2(7 are useful. We refer to this second group 
of three independent collaborations as Group B. 

What has made a direct comparison of the works of 
these papers impossible is that Group A used peak counts 
in maps of convergence, while Group B used peak counts 
in shear maps and a somewhat more realistic simulation 
of galaxy shape noise. For the first time in this paper, 
we use reduced shear maps and include realistic LSST 
measurement errors for the same set of simulations used 
by Group A, which allows a direct comparison of the 
works. However, it should be noted that our comparison 
is not complete, since we use only one smoothing scale 
and a different aperture mass filter compared to the work 
in Group A. 

We plot the A% 2 between different cosmological mod- 
els coming from the different SNR ranges in Figure [TUJ 
Neglecting correlations between individual SNR ranges 
and simply interpreting the area under the curves as the 
strength of distinction between the cosmological models, 
we conclude that peaks with SNR^ — 2a carry approx- 
imately 1/3 of all the information in the peak counts, 
peaks with SNR> 3rr approximately half, and peaks with 
S NR> 3.5cr also about 1/ 3. We can compare this result 
to IKratochvil et all (|2010l ). which found according to the 
third panel of their Figure 5 that low peaks (by our def- 
inition of 'low') were somewhat more important for cos- 
mology with convergence maps. However, we use a filter 



for aperture mass that emphasizes an NFW profile and 
so may de-emphasize smaller peaks, which would explain 
the discrepancy in our 

results^ 

In the literature, there have been claims that low (0- 
3(7 ) peaks do not carry any useful cosmological informa- 
tion both due to galaxy shape noise dominating the peak 
counts, and due to t he unknown i nfluen ce of systematic 
errors in this range. I Yang et"ail (|201lD has shown that 
the first of these issues is a misconception, and that real 
cosmological structure contributes significantly to peak 
counts at low SNR. We have shown here that, even in 
the presence of systematic errors from a realistic analy- 
sis pipeline, we can still extract that information. 




4 5 

SNR threshold 

Figure 10. The A% 2 (parameter sensitivity) from peak counts 
for our non-fiducial cosmological models, shown as a function of 
peak height in units of signal-to-noise ratio. Neglecting correlations 
between the different heights, positive peaks with SNR < 2a carry 
approximately 1/3 of the information content, those with SNR> 3cr 
carry one half, while peaks with SNR> 3.5a carry about another 
third. See Table [T] for details of cosmological model parameters. 



Inte restingly, an important result from Yan g et al.l 
(|2012|) is that peaks with SNR^ 1 — 3.5(7 are largely 
unbiased by baryon effects and therefore lend themselves 
particularly well for cosmological parameter estimations 
at the sub-percent accuracy level to which LSST aspires. 
Their result is valid within the restrictions of their study 
of varying the concentration parameter within NFW ha- 
los. To obtain certainty that this is universally the case 
for all possible contributions of baryonic physics, full hy- 
drodynamic simulations with different baryon prescrip- 
tions need to be run. We are in the process of investi- 
gating this for a future publication. 

We conclude that there is substantial information con- 
tent in low WL peaks, with 'low' defined as SNR~ — 3(7. 
Low peak counts should not be dismissed as purely due 
to shape noise, especially since these low peaks have been 
shown to be less sensitive to uncertain baryonic physics 
in simulations. 

6. ANALYSIS 

We present in this section the methodology for extract- 
ing cosmological information from the peak counts and 
power spectra we have obtained from aperture mass cal- 
culations. 



16 Also sec Pircs et al. (2012) which claims that convergence 
can contain complementary information to shear if one manages to 
extract it observationally. 
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6.1. Statistical Descriptors 

We generically refer to the different statistics one can 
obtain from a 2D WL map — e.g. power spectrum, 
peak counts, etc. — as statistical "descriptors", and de- 
note them by N. We can also combine them into a sin- 
gle vector, N{, where i indexes the peak height or the 
multipole for the power spectrum. Combining the data 
from several source redshifts makes the descriptor vector 
longer, but it is treated in the same way. We divide the 
range of peak height into 30 threshold bins. Similarly, 
we divide the power spectrum into 30 scale bins. 

To constrain cosmology, we are interested in the true 
ensemble averag43 (denoted henceforth by brackets ( )) 
and covariance of these descriptors as a function of cos- 
mological parameters (p = {O m , w, o~%}). These of course 
are not available to us, but can be estimated from the 
simulations. Averaging over the pseudo-independent 
map realizations within a given cosmology, we can es- 
timate the ensemble average by 

_ 1 R 

(Ni(p)) m Ni(p) = -Y,Ni(r,p), (4) 

r— 1 

where iVj(p) is the descriptor vector for one set of cos- 
mological parameters p, Ni(r,p) is the descriptor vector 
for a single realization and r runs over our R = 500 map 
realizations. We call this estimate the simulation mean. 
It differs from the true ensemble average both because 
of the limited number of realizations and also because 
of the limitations inherent in our simulations. In the 
absence of a fitting formula for the peak counts in the 
non-Gauss ian case (analogou s to the power spectrum for- 
mula from Smith et al. 2003) the simulation mean serves 
as our proxy for theoretically predicted peak counts P 5 ! 

Because of the computational expense of producing 
cosmological simulations, we can only produce a lim- 
ited number at the selected cosmologies given in Table [TJ 
We have calculated the simulation mean at these points 
and must extrapolate to other cosmologies not explic- 
itly simulated. Using finite differences between the sim- 
ulated cosmologies, we construct a first-order Taylor ex- 
pansion around our fiducial cosmology to estimate iVj(p) 
for other cosmologies: 

Ni(p) « Ni(p ) + ^ (p a - p 0a ), 

a Pa - POa 

Here, index a = 1,2,3 refers to an individual parame- 
ter, such as Q m ,w, or cr 8 . p( Q ) denotes the cosmological 
parameter vector of a simulated non-fiducial cosmology 
(where only the parameter p a has been varied) , and p 
denotes the parameter vector for the fiducial cosmology. 

The fraction in Eq. ((SJ is the finite difference deriva- 
tive. If the non- fiducial cosmology is chosen such that 

pi"' — po a is positive, we call it a "forward derivative" , 
if it is negative, we call it a "backward derivative". We 
compare the parameter constraints calculated from each 

17 Averaged over all possible universes with the same cosmolog- 
ical parameter values. 

18 While there are theoretical predictions for peak counts in 
Gaussian and non-Gaussian cases, the non-Gaussian predictions 
are not accurate enough to be useful for this purpose. 



derivative to assess the robustness of our results. 

Similarly to the simulation mean, we estimate the co- 
variance of the statistical descriptors from the simula- 
tions, Cov(Ni, Nj) s» dj, where 

1 R - _ 

Cij (p) = Yl w ^ p) - N * (p)i m ^ p) - ^ (p)i • 

r— 1 

(6) 

This covariance matrix contains contributions both from 
the sample variance of the true aperture mass signal and 
from the noise contributions. When the size of this co- 
variance matrix is large, inaccuracies in its estimate can 
become challenging, as we explore further below. 

6.2. Monte Carlo Probability Contours 

Each of our WL maps spans a 12 deg 2 field of view, yet 
we wish to obtain parameter contours for the full 20,000 
deg 2 LSST survey volume. We therefore employ boot- 
strapping to generate approximations to full-sky maps. 
In this procedure, we draw a map 20,000/12 ps 1667 
times from our 500 aperture mass maps, with replace- 
ment. The resulting 20,000 deg 2 map is not a true com- 
posite: we do not place the drawn maps edge to edge, 
but rather compute the descriptor values for each patch 
individually and then average over them to get their val- 
ues for the full-sky map. Details of this method, as well 
as its advantages for pa rameter estimation, will b e dis- 
cussed in this context in IKratochvil et aLl (|in prepl ). We 
create 10,000 such full-sky maps to obtain smooth pa- 
rameter contours in our Monte Carlo procedure. 

To estimate the cosmological parameter error contour 
from a set of WL maps from one cosmology, we use x 
minimization to fit for the best-fit cosmological parame- 
ters for each of the above bootstrapped full-sky maps. 
Thus our whole set of maps provides an ensemble of 
Monte Carlo realizations, and the distribution of those 
best-fit points can be used to draw probability contours 
at desired confidence levels. 

For realizations drawn from the fiducial cosmology p , 
X 2 is 

X 2 (r,p) =J2 A^(r,p)[Cov- 1 (p )] lJ A7V,(r,p) (7) 
where 

AN i (r,p)=N i (r,p )-{N i (p)). (8) 

For each Monte Carlo realization, we minimize x 2 with 
respect to p using a simulate d anneali ng algorithm. As 
outlined in IKratochvil et al.l ([in prepl ). the covariance 
matrix is computed from the auxiliary map set and in- 
verted with singular value decomposition, discarding any 
problematic eigenvectors. The simulation mean for the 
fiducial cosmology is computed from the primary map 
set, whereas the finite difference derivatives are com- 
puted from the auxiliary map set. The maps for which 
best-fit parameters are computed come from the primary 
map set. 

7. COSMOLOGICAL CONSTRAINTS 

In this section we present the results obtained by ap- 
plying the methods described in Section [6] to the peak 
counts described in Section [5] in order to constrain the 
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cosmological parameters f2 m , as and w. We also calcu- 
late the power spectrum of the aperture mass maps, and 
use that to predict co nstraints on cosm ological parame- 
ters. Previous work (jYang et all 1201 il ) has shown that 
peak counts in convergence maps contain additional in- 
formation beyond the power spectrum, and can tighten 
constraints on cosmological parameters, as well as break 
degeneracies in constraints from the power spectrum 
alone. We compare here the constraints on cosmologi- 
cal parameters obtained from the power spectrum of the 
aperture mass maps, traced by galaxies with shape noise 
included, and the constraints from peak counts, traced 
by the same galaxies. 

To determine the sensitivity of peak counts to cosmo- 
logical parameters, one can use the backward, forward, or 
symmetric derivative in the Taylor expansion in Eq. ([5])- 
It is important to check that all derivatives give the same 
contours, and as expected we find very small shifts in the 
contour size and shape between the different derivative 
methods consistent with the statistical limitations of the 
cosmological simulations used. The contours we present 
in this section are 68% error contours, corresponding to 
data obtained by a full LSST ten-year survey using only 
good quality r-band data. 

Figure [11] shows the contours for the peak counts above 
SNR thresholds for the backward derivative. The con- 
tours in dashed curves show the predicted constraints 
for measurements made with shape noise only, and solid 
curves show the constraints for both intrinsic shape noise 
and the realistic measurement errors described in previ- 
ous sections. The predicted constraints with measure- 
ment error appear at first inspection to be smaller than 
with shape noise alone. In fact this is due to the ac- 
curacy with which we can derive the contours from the 
available set of cosmological simulations. A set of simu- 
lations covering a larger range of cosmological parameter 
space would yield smoother contours, but at present this 
is computationally unfeasible. If we look at the contours 
calculated for the forward derivative shown in Figure 1121 
we see that in this case the contours with measurement 
error are roughly the same size as those for shape noise 
alone. The constraints with and without measurement 
error are indistinguishable within the statistical accuracy 
of the contours. 

This indicates that measurement errors will have a rel- 
atively small impact on the accuracy of cosmological con- 
straints with shear peak statistics with LSST. We note 
that we have neglected spatial correlations in measure- 
ment errors, which we expect to have a small impact on 
the number counts of peaks but are known to have a 
more significant impact on measurements of the power 
spectrum. In comparing constraints obtained with peak 
counts and our measurement of the power spectrum, we 
therefore restrict the comparison to the case with shape 
noise only. This is shown in Figu r e IT51 As seen in p revi- 
ous work (jKratochvil et all 1201% lYang et all 1201 II ). the 
predicted constraints obtained with peak counts are bet- 
ter than using the power spectrum alone, and a small 
improvement is found when combining the two measure- 
ments. 

It should be noted that our finding that measurement 
errors have a small impact on constraints from peak 
counts, is only valid under the assumptions that we have 
made in our analysis framework. For example, there 



may be a larger error associated with the KSB algorithm 
than we find in this work, because our simulated galaxies 
are modeled with Sersic profiles rather than real galaxy 
shapes. Since the ImSim input catalogues are anchored 
to real data, this is a limitation from the current sur- 
vey data, and will be improved with future observations. 
The accuracy of our constraints is also limited by the ac- 
curacy of our N-body si mulations, which i s disc ussed in 
detail in Section V.F of Kratochvil et all (|2012l) . 

The predicted constraints in this work are comparable 
to those found in other analyses which use the same cos- 
mologi cal simulations, but different measurem ent tech- 
niques (|Kratochvil et alJl2012tlYang et al.|[201ll) . For ex- 
ample, Figure 12 in lYang et al. I (|2011l) shows constraints 
on cosmological parameters from peak counts from con- 
vergence in a single redshift plane, combined with the 
power spectrum and scaled to an LSST-size survey, that 
are very similar to our constraints[3 The comparison 
is not direct, however, due to the differences in analysis 
methodology. In particular, most of the previous work 
has looked at peaks in maps of convergence, whereas we 
study reduced shear p eak counts. The s imilarity of our 
constraints to those in Ya ng et all (|2011l ) imply that the 
difference between peak counts in convergence and re- 
duced shear is small, but we are unaware of any work 
directly comparing the two methodologies. 

It is even harder to make a comparison with other work 
that uses an entirely different set of cosmological simu- 
lations. In that case, the cosmological parameters varied 
in the simulations, as well as the details of the N-body 
algorithms, make it almost impossible to make meaning- 
ful comparisons between results. The best we can say is 
that the amount of information added by a measurement 
of peak counts, compared to the power sp ectrum alone, 
is consistent with results from other work (|Marian et all 
ll^lDietrich fc Hartlaijll20Toh . 

It should be noted that this work does not use multi- 
ple smoo thing scales. Using Mink owski functionals as an 
example, iKratochvil et all (|2012l) have shown that com- 
bining smoothing scales is important to extract the max- 
imum amount of information from weak lensing maps 
with a non-Gaussian descriptor, and lMarian et all()2012[ ) 
have explicitly demonstrated this for peak counts. It re- 
mains to be studied if the power spectrum contributes 
any additional information to peak counts when a com- 
bination of smoothing scales is used for peaks, or if peaks 
manage to extract all informatio n when enough smooth - 
ing scales are used. The work in IKratochvil et all (|2012f ) 
noted that combining smoothing scales does not improve 
the constraints of the power spectrum, because the small- 
est smoothing scale contains already all of the informa- 
tion contained in the power spectrum. This is true for the 
Gaussian filter used on convergence in that paper. How- 
ever, when evaluating aperture mass on reduced shear, 
one arrives at a compensated filter on convergence. Com- 
pensated filters suppress all modes much longer than the 
size of the filter, in addition to small-scale modes. Thus, 
for the filter used in this paper — and for using shear in 
general — combining different smoothing scales could be 
an asset also for the power spectrum. 

19 It should be noted that the analysis in Yang et al. ( 2011) used 
only 15 galaxies per arcminute 2 for one redshift plane, whereas we 
use 30 galaxies per arcmin 2 over the full redshift range. 
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Figure XI. 68% error contours on the cosmological parameters r2 m , erg and w for peak counts, using the backward derivative in the Taylor 
expansion. The dashed curves show the constraints for measurements including shape noise only, and the solid curves for both shape noise 
and measurement error. 
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Figure 12. 68% error contours on the cosmological parameters f2 m , erg and w for peak counts, using the forward derivative in the Taylor 
expansion. The dashed curves show the constraints for measurements including shape noise only, and the solid curves for both shape noise 
and measurement error. 
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Figure 13. 68% error contours on the cosmological parameters fi m , erg and w for peak counts (blue) and power spectrum (red) and the 
combination of the two (green) , using the backward derivative in the Taylor expansion for measurements including shape noise only 



8. SUMMARY 

We have produced the first framework for including 
realistic galaxies and measurement errors in predictions 
of shear peak counts, using information from the LSST 
Image Simulator, ImSim. Galaxies are drawn from real- 
istic distributions, based on observational data, in red- 
shift, size, magnitude and ellipticity. We use information 
from ImSim to assign uncertainties to the galaxy shape 
measurements based on these properties, using the KSB 
shape measurement algorithm. We use these realistic 
galaxies to trace the reduced shear maps produced from 
ray-traced cosmological N-body simulations, and distort 



the galaxy shapes appropriately according to the shear 
parameters interpolated in three dimensions. The aper- 
ture mass and signal-to-noise ratio is calculated for the 
resulting simulated catalogues, using an implementation 
of the aperture mass statistic on the GPU. 

We count peaks above SNR thresholds, and use the 
resulting peak count distributions to predict constraints 
on cosmological parameters with LSST. We also calcu- 
late the aperture mass for an idealized case where we 
know the reduced shear perfectly, with no uncertainty 
from galaxy shape noise or measurement error. Compar- 
ing the two cases, we find that the majority of the dis- 
criminating power for the ideal case is in the high SNR 
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peaks, whereas for the realistic measurements the power 
comes from the full range of peak counts. This confirms 
for the case of peak counts in reduced shear maps what 
has alrea dy been seen in p eak c o unts for convergen ce 
maps by iKratochvil et all (|2010f ): lYang et al.1 (|2Qlll )— 
that low and medium significance peaks with SNR< 3.5<r 
in reduced shear maps contribute most of the cosmolog- 
ical constraining power. 

We calculate the 68% confidence contours for the re- 
alistic and noiseless peak counts, and find that there is 
a significant degeneracy in Sl m and ag, and smaller de- 
generacies in the other planes. Even in the presence of 
noise, there is substantial information in peak counts be- 
yond that which can be extracted from the power spec- 
trum alone - the contours from peak counts are approxi- 
mately 25% smaller compared to contours from the power 
spectrum alone. Combining the two measurements (from 
peak counts and the power spectrum), we calculate ex- 
pected constraints on fi m , as and w for a ten- year LSST 
survey of 0.257 < Q m < 0.263, 0.792 < cr 8 < 0.804 and 
0.98 < w < 1.02. Note that these constraints are given 
for illustration, and are not intended as expected final 
constraints for LSST. LSST analyses will use multiple 
redshifts bins for power spectrum and peaks and include 
measurements of BAO, supernovae etc. to achieve sub- 
stantially better constraints. 

We have shown that reduced shear peak counts are 
a useful probe of cosmology, and that the presence of 
realistic instrument noise and measurement uncertainties 
have very little impact on the power of the cosmological 
constraints. This work is the first step in a series of 
analyses needed to develop the analysis of shear peak 
counts to the level of sophistication currently enjoyed 
by the study of shear correlation functions. Future work 
should consider the impact of photometric errors, masked 
areas, and galaxy clustering around areas of high shear 
on peak counting, and should start to consider ways to 
potentially mitigate these effects. 

We would like to thank Brandon Calabro for useful 
discussions on shear maps and aperture mass statistics. 
J.M.K. and K.M.H. receive support from NASA's Jet 
Propulsion Laboratory subcontract 1363745. This re- 
search utilized resources at the New York Center for 
Computational Sciences, a cooperative effort between 
Brookhaven National Laboratory and Stony Brook Uni- 
versity, supported in part by the State of New York. This 
work is supported in part by the U.S. Department of En- 
ergy under Contract No. DE-AC02-98CH10886. 

REFERENCES 

Allen, S., Evrard, A. & Mantz, A., 2011, ARAA 49, 409. 

Albrecht A. et al., 2006, |arXiv:astro-ph/060959ll 

Bard, D., Bellis, M., Allen, M. et al., 201!2,|arXiv:1208.3658vl 
[astro-ph.IM]. 

Bartelmann, M., & Schneid er, P., 2001, Ph ysics Reports, 340 291. 
Behroozi, P. et al., 2011. larXiv:1110.4370fr 2 [astro-ph.CO] . 
Bertin, E., Arnouts, S., 1996, AAS 117, 393. 
Bruzual, G. & Chariot, S., 2003, Mon. Not. Roy. Astron. Soc. 
386, 887. 

Coil, A. et al., 2004, Astrophys. J. 609, 525. 
Chang, C, et al., 2012, arXiv:1206.1378v2 [astro-ph.CO]. 
Chang, C, et al., 2012, arXiv:1206.1383v2 [astro-ph.IM]. 
Connolly, A., et al., 2010, Simulating the LSST System, Proc. 
SPIE 7738, 53. 



De Lucia, G. et al., 2006, Mon. Not. Roy. Astron. Soc. 366, 499. 
Dietrich, J. P. & Hartlap, J., 2010, Mon. Not. Roy. Astron. Soc. 
402, 1049. 

Guimara.es, A., 2002, Mon. Not. Roy. Astron. Soc. 337, 631 . 
Hamana, T. & Mellier, Y., 2001, Mon. Not. Roy. Astron. Soc. 
327, 169. 

Hamana, T., Takada, M. & Yoshida, N., 2004, Mon. Not. Roy. 

Astron. Soc. 350, 893. 
Hetterscheidt, M., Erben, T. & Schneider, P., 2005, A & A 442, 

43. 

Hockney, R. & Eastwood, J., 1998, "Computer Simulation Using 

Particles," Adam Hilger, Bristol. 
Hoekstra, H. & Jain, B., 2008, Annual Review of Nuclear and 

Particle Science 58, 99. 
Jain, B. & Van Waerbeke, L., 2000, Astrophys. J. Lett. 530, LI. 
Jain, B., Seljak, U. k. White, S., 2000, Astrophys. J. 530, 547. 
Jain, B., Seljak, U. & W hite, S., 2000, Ast rophys. J. 530, 547. 
Joachimi, B. ct al., 2012. larXivTi 203.6833vl [astro-ph.CO]. 
Kaiser, N., Squires, G. & Broadhurst, T. 1995, Astrohphys. J. 

449, 460. 

Kaiser, N., 1995, Astrophys. J. 439, LI. 

Kitching, T. et al., 2012, Mon. Not. Roy. Astron. Soc. 423, 4. 
Kratochvil, J. M., Haiman, A. & May, M., 2010, Phys. Rev. D 
81, 043519. 

Kratochvil, J. M., Wang, S., Lim, E. A. et al., 2012, Phys. Rev. D 
85, 103513. 

Kratochvil, J. M., "The Inspector Gadget Weak Lensing 

Simulation Pipeline: I. The Codes" , in prep. 
Kratochvil, J. M., May, M., Haiman, Z. et al., "The Inspector 

Gadget Weak Lensing Simulation Pipeline: II. The First 

Dataset" , in prep. 
Leauthaud, A. et al., 2007, Astrophys. J. S 172, 1. 
Lewis, A., Challinor, A., & Lasenby, A., 2000, Astrophys. J. 538, 

473. 

Li, M., 2011, A & A 35, 349. 

Lombardi, M., & Bertin,, G., A & A 335, 1 (1998). 
Marian, L., Smith, R. E., & Bernstein, G. M., 2009, Astrophys. J. 
698, L33. 

Marian, L., Smith, R. E., & Bernstein, G. M., 2010, Astrophys. J. 
709, 286. 

Marian, L., Hilbert, S., Smith, R. E. et al, 2011, Astrophys. J. 
728, L13. 

Marian, L., Smith, R. E., Hilbert, S. et al, 2012, Mon. Not. Roy. 

Astron. Soc. 423, 1711. 
Maturi, M., Angrick, C, Bartelmann, M. et al, 2010, A & A 519, 

A23. 

Maturi, M., Fedeli, C, & Moscardini, L., 2011, Mon. Not. Roy. 

Astron. Soc. 416, 2527. 
Munshi, D., Valageas, P., van Waerbeke, L. et al., 2008, 

Phys. Rep. 462, 67. 
Munshi, D.„ van Waerbeke, L., Smidt, J. et al., 2011, Mon. Not. 

Roy. Astron. Soc. 419, 536. 
Navarro, J., Frenk, C. & White, S., Astrophys. J. 486, 493 (1996). 
Peterson, J. et al., 2009, LSST Science Book, Chapter 3.3, 

arXiv:0912.0201 , http://www.lsst.org/lsst/scibook 
P into, P. et al., 20097LSST Science Book, Chapter 3 . 1 , 

arXiv:0912.0201 . http://www.lsst.org/lsst/scibook 
Pires, S., Starck, J-L., Amara, A. et al., 2009, A&A505, 969. 
Pires, S., Leonard, A., & Starck, J-L., 2012, Mon. Not. Roy. 

Astron. Soc. 423, 983. 
Sato, J., Takada, M., Jing, Y. P., & Futamase, T., 2001, 

Astrophys. J. 551, L5. 
Sato, M., Hamana, T., Takahashi, R., et al., 2009, Astrophys. J. 

701, 945. 

Schirmer, M., et al., 2007, A+A 462, 875. 

Schneider, P., Ehlers, J., & Falco, C. C, 1992, "Gravitational 

Lenses" , Springer- Verlag, New York. 
Schneider. P., 2005, |arXiv:astro- ph/0509252 
Schrabback, T, et al., 2010, AfeAS16, £5X 
Seitz, C, Schneider, P., 1998, |arXiv:astro- ph/9802051 vl. 
E. Semboloni et al., 2011, Mon. Mot. Roy. Astron. Soc. 410, 143. 
Smith R.E., Peacock, J. A., Jenkins, A. et al., 2003, Mon. Not. 

Roy. Astron. Soc. 341, 1311. 
Springel, V., 2005, Mon. Not. Roy. Astron. Soc. 364, 1105. 



14 



Bard et al. 



Turner, E. L., Ostriker, J. P., Gott, J. R., 1984, Astrophys. J. 
284, 1; Webster, R. L., Hewett, P. E., Harding, M. E., et al., 
1988, Nature 336, 358; Fugmann, W., 1988, Astron. & 
Astrophys. 204, 73; Narayan, R., 1989, Astrophys. J. Lett. 339, 
53; Schneider, P., 1989, Astron. & Astrophys. 221, 221. 

VanderPlas, J., et al. 2012, Astrophys. J. 744, 180. 

Wambsganss, K., Cen, R., & Ostriker, J. P., 1998, Astrophys. J. 
494, 29. 

Wang, S., Haiman, Z., & May, M., 2009, Astrophys. J. 691, 547. 



Wittman, D. et al., 2000, Nature 405, 143. 

Witt man, D. et al. , 2009, LSST Science Book, Chapt er 14, 

arXiv:0912.0201 http://www.lsst.org/lsst/scibook 
Yang, X., Kratochvil, J.M., Wang, B. et al., 2011, Fhys. Rev. D 

84, 043529. 

Y ang, X., Kratochv il, J. M., Huffenberger, K. M. et al., 2012, 
arXiv:1210.0608 astro-ph.CO]. 



